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Abstract 

We consider a one-dimensional continuous time random walk (CTRW) on a fixed time interval T 
where at each time step the walker waits a random time r, before performing a jump drawn from 
a symmetric continuous probability distribution function (PDF) /(r/), of Levy index 0 < p < 2. 
Our study includes the case where the waiting time PDF 'k(r) has a power law tail, 'l'(r) oc , 

with 0 < 7 < 1, such that the average time between two consecutive jumps is infinite. The random 
motion is sub-diffusive if 7 < p/2 (and super-diffusive if 7 > p/2). We investigate the joint PDF 
of the gap g between the first two highest positions of the CTRW and the time t separating these 
two maxima. We show that this PDF reaches a stationary limiting joint distribution p{g, t) in 
the limit of long CTRW, T —)• 00 . Our exact analytical results show a very rich behavior of this 
joint PDF in the ( 7 , p) plane, which we study in great detail. Our main results are verified by 
numerical simulations. This work provides a non trivial extension to CTRWs of the recent study 
in the discrete time setting by Majumdar et al. {J. Stat. Mech. P09013, 2014). 
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I. INTRODUCTION AND SUMMARY OF MAIN RESULTS 


Recently, there has been a resurgence of interest in extreme value questions for random 
walks and its variants, particularly in the physics literature [1-6]. This is motivated to a 
large extent by the observation that extreme value statistics (EVS) plays a crucial role in 
the physics of complex and disordered systems [7, 8 ]. Although the EVS of independent 
and identically distributed (i.i.d.) random variables is very well understood, much less is 
known for strongly correlated variables, which is the case frequently encountered in physics. 
From that point of view, random walks offer a very useful framework where the effects of 
correlations on EVS can be characterized analytically. 

The statistics of the maximal displacement of a random walk (RW) after a given number of 
steps n is now quite well understood, thanks for instance to the Pollaczeck-Spitzer formula [9, 
10], from which precise asymptotic estimates in the large n limit can be obtained [1]. The 
full order statistics, i.e. the statistics of the maximum for any integer k, was investigated 
in [11] in the case of ordinary RWs with a narrow jump distribution (characterized by a Levy 
exponent p = 2 ). Of particular interest is the gap between two consecutive maxima, which 
makes it possible to characterize the “crowding” of the walker positions near their maximum. 
Along this line we recently obtained a complete description of the joint distribution of the 
gap and time interval between the hrst two maxima of RWs not only for narrow jump 
distributions (p = 2), but also for Levy flights characterized by a Levy exponent 0 < p < 2 
[4, 5]. Since the characteristic displacement of the walker after n time steps grows like 
the results of [4, 5] are of interest for diffusive (/i = 2) and super-diffusive (0 < p < 2) 
processes only, leaving out the important case of sub-diffusive processes. Extending the 
analysis of [4, 5] to sub-diffusive RWs has been one of the motivations for this work. 

A common model to describe sub-diffusive processes is the so called continuous time 
random walk (CTRW). In this model, the walker performs a usual random walk where 
two consecutive jumps are separated by a certain trapping time r. The trapping times are 
i.i.d. random variables drawn from a common distribution T(r) which has a power law tail 
^/(r) oc with 0 < 7 < 1. This model, initially suggested by Scher and Montroll in the 

setting of non-Gaussian transport in disordered electronic systems [12], has been widely used 
to describe anomalous sub-diffusive dynamics in various complex systems phenomenologi¬ 
cally [13, 14]. Indeed, for 0 < 7 < 1 the average trapping time between two successive jumps 
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is infinite and, consequently, CTRWs are characterized by an anomalous growth of displace¬ 
ment oc which is indeed sub-diffusive if 7 < /r/2 (and super-diffusive for 7 > /i/2). The 
case of a hnite average trapping time between two successive jumps corresponds to standard, 
non sub-diffusive, RWs. The fluctuations of the maximal displacement of CTRWs for /i = 2 
and arbitrary 7 < 1 was studied in Ref. [ 6 ] using real space renormalization techniques 
(see also [15] for a Feynman-Kac approach). A precise estimate of the expected value of the 
maximum of the CTRW was obtained in Ref. [16] by combining the results of Ref. [1] with 
the so called ’subordination property’. However, nothing is known about the order statistics, 
beyond the first maximum, for such CTRWs. In this paper, we provide a complete analytic 
description of the distribution of the gap and the time interval between the hrst two maxima 
of such CTRWs for any values of the parameters 0 < /r < 2 and 0 < 7 < 1, as well as for a 
hnite average trapping time between two successive jumps. 

We consider a CTRW starting at the origin, x(t = 0) = Xq = 0, and evolving according 
to 

Xi = Xi_i + r]i, ( 1 ) 

where Xi denotes the walker position between the Tth and the {i + l)-th jumps. The jumps 
r]i are i.i.d. random variables distributed following a symmetric, bounded and piecewise 
continuous distribution /(?/) the Fourier transform of which, f{k) = has 

the small k behavior 

f[k) = l-\ak\^ + o{\k\n, (2) 

where 0 < /r < 2 is the Levy index and a > 0 is the characteristic length scale of the jumps. 
For n = 2, the variance of the jump distribution rff{rj)drj is well dehned and 

a = (j/\/2. On the other hand, for 0 < /i < 2, /(?/) does not possess a well dehned second 
moment because of its heavy tails, /(?/) oc {rj —)■ 00 ), and the RW (1) is a Levy hight 

of index /r. The time intervals between the {i — l)-th and the f-th jumps (r* > 0) are 
i.i.d. random variables, independent of rji, with PDF 'L(r) the Laplace transform of which, 
'F(( 3 ') = ^'(t) exp(—gr) dr, has the small q behavior 

'^( 9 ) = 1 - (R9)^ + o(|gr), (3) 

where 0 < 7 < 1 and Tc > 0 is the characteristic time scale of the jumps. For 7 = 1 , the 
mean time between two successive jumps, (r) = r\k(r) dr < -|-oo, exists and Tc = (r). 
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For 0 < 7 < 1, the mean time between two successive jumps does not exist. As a function 
of time, the walker position is thus given by 


dx(t) 

dt 




(4) 


i=l \ j=l / 

with a;(0) = 0 and n the total number of jumps in the walk. The solution to (4) is readily 
found to be given by F?' < ^ < Si=i F?) ~ with Xi given by ( 1 ). 

The goal of the present paper is to provide an exhaustive discussion of the different 
behaviors of the joint PDF of the gap g and time interval t between the first two maxima 
of the walk that may arise in the large n limit, depending on the large argument behavior 
of f{rf) and 4/(t). Before entering the details of the calculations, it is useful to summarize 
our main results. We first show that the joint PDF Pn{g,t) has a well dehned limiting PDF 
as n —)■ cx): 


\im pn{g,t) = p{g,t), (5) 

n^oo 

where the Laplace transform of p{g,l) with respect to (w.r.t.) t is given in Eqs. (20) and 
(17). Then, we perform a detailed analysis of p{g,t) in the plane {g,t) for p and 7 in the 
whole ranges 0 < /i < 2 and 0 < 7 < 1 , and for the three different main classes identihed in 
[5]: (i) slow, (ii) exponentially, and (hi) fast decreasing f{rj) at large 7 . 

One of the most remarkable result of this study may be the existence of a scaling form 
for the joint PDF p{g, t) in the case of a fast decreasing jump distribution (/i = 2), such that 
fiv + x) f (p) exp{—cp^ x) for 7 —)■ cx) [see Eq. ( 68 ) below for a more precise dehnition], 
with 0 < 7 < 1. This class includes the case of a Gaussian jump distribution for which 
(5 = 1. This is a new result, characteristic of CTRWs, without any counterpart in the 
discrete time RW setting considered in [5]. More specihcally, we show that there is a scaling 
regime g, 1 with hxed in which p{g,t) takes the scaling form 

(9 ^+00 and i ^+00), (6) 

with the asymptotic behaviors 


K{y)^ 


Pi^-2(i+i/7)^ (y^O) 

{y+ 00 ) 


(7) 


where the amplitudes Vi and Vu are given in Eqs. (80) and (81), respectively. Physically, 
the switch from the hrst to the second behavior (7) around ~ acg^ corresponds to 
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the cross-over from a ‘concentration’ - or ‘one-step’ ~ regime (for < acg^) where 

the walker get stuck for a long time t at the second maximum and then jumps directly to 
the hrst maximum, to a ‘many-steps’ regime (for (f/rc)'’'^^ > acg^) where she/he travels a 
long walk of total duration t (with many steps) between the second and the hrst maxima. 
For 7 = 1, there is no scaling form and ( 6 ) is replaced with the uniform expression (84). 

An other important result is how the scaling form derived in [5] for a discrete time Levy 
hight (see Sec. 4.1 and Appendix D in [5]) is affected when one switches to a CTRW. 
Levy hights have type (i) jump distributions since, for 0 < /i < 2, /(?]) necessarily has an 
algebraic, slow decreasing, tail. In this case we show that in the scaling regime g, t ^ 1 
with hxed Vg~^, the joint PDF p{g,t) takes the scaling form 

t) -( - ) F /,7 ( —^ ) {g +00 and t +cx)), (8) 

arc \gj J 

with the asymptotic behaviors 

Cj I _ int(7) < ^ < 2, 

< Cuy-^-F^\ny, p = l, {y ^+oo), (9) 

0</i<l, 

\ 

and 

F,^,{y)^CjyyF^-F\ (y^O), ( 10 ) 

where Int( 7 ) in the hrst line of (9) denotes the integer part of 7 and the last two lines are 
for 0 < 7 < 1 only. It can be checked that for 7 = 1, Eqs. ( 8 ), (9), and (10) reduce to Eqs. 
(87), ( 88 ), and (90) of [5] in which the number of jumps between the hrst two maxima, /, 
is merely replaced with t/rc, as expected from simple law of large number arguments giving 
I —)■ tjTc as f —)■ -1-cx). Our results provide thus a non trivial extension of the ones in [5] to 
the case 0 < 7 < 1 in which the mean time between two successive jumps does not exist 
and law of large number arguments cannot be used. 

Finally, by integration of the joint PDF p{g,t) over g, one hnds that the marginal dis¬ 
tribution ptime(^) displays a power law tail with logarithmic corrections and an exponent 
depending only on p, and 7 . Namely, for 7 = 1 one hnds 

f ^I(^/re)-l-l/^ l</i<2, 

Ptime(f) ~ — X < An ln(t/rc) h = + 00 )> (11) 

x4iii (t/rc)“^ 0 << 1, 
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where the amplitudes Aj, An and are given in Eq. (91), and for 0 < 7 < 1 one has 

Ptimeit) ~ — X < ln^(t/rc) /i = 1 , (f +00), (12) 

^ni ln(t/rc) {t/Tc)~^~'^, 0 < /i < 1, 

\ 

where the amplitudes A!i, A{i and are given in Eq. (97). Again, it can be checked that 
Eq. (11) is nothing but Eq. (112) in [5] with / = t/rc, as expected (c/. the end of the 
preceding paragraph). Note that for /i = 2 and 7 = 1, one has ptime(^) oc whatever the 
distributions f{r]) and '^{t) possessing a second and a first moment, respectively. The third 
lines of ( 11 ) and ( 12 ) reveal an interesting freezing phenomenon, as a function of n, of the 
large t behavior of ptime(t) as fx decreases past the value fXc = ^- It follows in particular from 
(11) and (12) that the hrst moment of ptime(t) is never dehned. This means that, although 
the typical size of t is 0(1), its average diverges with the total duration, T, of the random 
walk. For 0 < 7 < 1 one can estimate from (11) and (12) that {\t\) ~ for 1 < /x < 2, 

while {\t\) ~ T^~^\ogT for 0 < /r < 1 and {\t\) ~ T^“'^(lnT)^ for n = 1. Table I summarizes 
the functional dependence of the different asymptotic behaviors of Ptime(^) at large t in the 
(/i, 7 ) plane. 


Ptime(t) OC 

1 < p < 2 

p = 1 

0 < p < 1 

7 = 1 

t-i-iA 

Int 

t-2 

0 < 7 < 1 


^-i- 7 (ln ^)2 

1 

7 


TABLE I: Functional dependence of the different asymptotic behaviors of ptime(t) at large t in the 
(p, 7 ) plane. 

Note that the marginal distribution of the gap between the first two maxima, Pga.p{g), 
obtained by integration of the joint PDF p{g,t) over t, does not depend on T(r) (hence on 
7 ) and is the same as the one already studied in Sec. 3 of Ref. [5] in the discrete time RW 
setting. 

The outline of the paper is as follows. In Section II we give the Laplace transform of p{g, t) 
with respect to t for a free-end random walk and a random bridge in which the walker is 
conditioned to return to a; = 0 at the end of the walk. Section III deals with the asymptotic 
behavior of p{g, t) for large t at hxed g, from which we derive, in Section IV, the asymptotic 
behavior of p{g, t) for a Levy flight at large t then large g. The asymptotic behavior of p{g, t) 
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for large g is investigated in Section V for three classes of jump distributions encompassing 
a wide range of jumps of practical interest. It is shown in Section V C that when the jump 
distribution has a fast decreasing tail, the behavior of p{g, 1) when both g and I are large takes 
on a scaling form capturing the existence of two complementary regimes, one dominated by 
walks with only one step between the hrst two maxima, and the other dominated by walks 
with many steps between the hrst two maxima. In Section VI we give the scaling form of the 
asymptotic behavior of p{g, t) for a Levy hight at large t and g. In Section VII we use the 
results obtained in Sections III and VI to determine the large t behavior of Ptime(^) and show 
its freezing when the random walk is a Levy hight with index 0 < /x < 1. Finally, Section 
VIII is devoted to the comparison of our analytical results with numerical simulations. 


II. LAPLACE TRANSFORM OF p{g, t) WITH RESPECT TO t FOR A FREE-END 
RANDOM WALK AND A RANDOM BRIDGE 


Let Pn{g, t) denote the joint PDF of the gap, g, and the time interval, t, between the hrst 
two maxima of a continuous time random walk after n steps. Since Pn{g, —t) = Pn{g,t) we 
will restrict ourselves to t > 0 without loss of generality. Let / denote the number of jumps 
between the hrst to maxima. For t > 0 one has 


Pn{g,t) = ^Pn{g,l,t)= ^ p{t\g,l)pn{g,l) = ^ p{t\l)Pn{g,l), 

0<l<n 0<l<n 0<l<n 

where p{t\l) has the Laplace transform 


r*+oo 


p{q\1) = / p(^l0 exp(-gt) dt = 


JJ (exp(-gri)) = 'F(g)^ 

i=l 


/*+oo / ^ 

/ 


Ti) ) exp(—gt) dt 


Taking the Laplace transform of (13) with respect to t and using (14) one gets 


^+00 


Pn{g,q) = / Pn{g,t) exp{-qt)dt= ^ 'L(g)'pn(^,/) = Pn(^,'L(g)), 


0 </<n 


(13) 


(14) 


(15) 


where Pn{g, s) = Ylio<i<nPn{9) 0 is the generating function of Pn(fl', 0 with respect to / > 0. 
The large n limit of (15) can then be readily obtained from the main result of [5] according 
to which p(g, s) = hm„^+ooPn(5', s) exists and is given by 


P{9,s) = /i(g,s)/2(g). 


(16) 
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with 


h{.g, s) = s uiy, s)f{g + y) dy, 
hid) = Jo°^ h{x, l)f{g + x) dx, 

where the functions u{x, s) and h{x, s) are the inverse Laplace transforms of 4>{X, s) and 
0(A,s)/A, respectively: 


Io°^'^{x,s)e ^^dx = s), 

h{x, s)e~^^dx = 0(A, s)/A, 


(18) 


with 


(j){X, s) = exp- 


A f+°^\n[l-sf{k)] 


TT 


A;2 + A2 


dk 


Injecting this resnlt into the large n limit of (15) yields 


(19) 


Pla.tl) = lim p..(j.9) =p(9. *( 9 )) = A(j,'i'(<l))/2(9). ( 20 ) 

n^+oo 

Thus, Pn{g,t) converges to a limiting distribntion p{g,t) as n —)■ +cx) whose Laplace trans¬ 
form with respect to t, p{g,q) = p{g,t) exp{—qt) dt, is given by ( 20 ). 

It was proved in [5] that the expression (16) of p{g, s) actnally holds for both free-end 
random walks and random bridges in which the walker is conditioned to retnrn to a; = 0 
at the end of the walk. Since, for a given T(r), p{g,q) is entirely determined by p{g,s) 
through ( 20 ), it follows immediately that the expression of p{g,q), hence the one of the 
limiting joint PDF p{g, t), is exactly the same for free-end walks and bridges. Therefore, all 
the results obtained in the following apply to free-end random walks as well as to random 
bridges withont any modihcation. 


III. ASYMPTOTIC BEHAVIOR OF p{g, t) FOR LARGE t AT FIXED g 

Sections III to VI are devoted to the asymptotic behavior of the joint PDF p{g, t) when 
either g or t (or both) is large. From the results obtained for the large t behavior of 
p{g, t) it will then be possible to derive the large t behavior of the marginal distribntion 
Ptime(^) = p{g,t) dg, which will be the snbject of Sec. VII. 

First, we consider the limit t —)■ -|-cx) at hxed g. As we will see, the behavior of p{g,t) in 
this limit depends on the valne of 7 , leading to new, non trivial, asymptotic expressions if 
0 < 7 < 1 . 
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A. 7 = 1 


In this case the mean time between two successive jumps, (r) = Tc, exists and by a simple 
law of large number argument, the large t behavior of p{g, t) can be expected to be given by 
T~^p{g,l = t/Tc), where I = tjTc is large. Let us check this result explicitly by comparing 
the large t behavior of p{g, t) and the large I behavior of p{g, /). 

The behavior of p{g, t) at large t is determined by the one of p{g, q) at small q. Thus, 
from (20) and (3) with 7 = 1, one has 


P{9,t) = Pia,(q)) exp{qt)dq 

~ ^ / P{9, 1 - Tcq) exp(gt) dq {t +oo) 


1 


2m T, 


Pi9 :1 - ?) exp ( — ) dq, 


C J C 


( 21 ) 


where we have made the change of variable q —)■ q/Tc- Similarly, for large I the behavior 
of p{g, 1) is determined by the one of p{g, s) near s = 1. Making the change of variable 
s = exp(—g) one gets 




P{9,s) 


ds 


|•+i'K+0 

= 7^ P(^,exp(-g)) exp(g/)dg 

' —-iTT+O 


2m 


^ / P{9, 1 - g) exp(g/) dq (/ ^ +cx)). 


( 22 ) 


Comparing Eqs. (21) and (22), one obtains the expected expression 


P(^,f)- p{9,l=-] (t^+oo). 


Tr 


Finally, Eq. (65) of [5] together with Eq. (23) above yield, for all 0 < p < 2, 

, , r(l + l//i) Ugf 

Pi9: t) ~ ^^ u, ^ +'^)- 

naTc \t/Tc) ^ 


(23) 


(24) 


B. 0 < 7 < 1 

If 0 < 7 < 1 the mean time between two successive jumps, (r) = +cx), does not exist and 
the large t behavior of p{9,t) is no longer related to p{9,l) by the simple expression (23). 
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Following the same line as for Eq. (21) with 0 < 7 < 1, one now has 


P{9,t) = ^ / Pia,(q)) exp{qt)dq 

~ ^ / p{9 ,1 - {rcqy) exp{qt) dq {t + 00 ) 


( 25 ) 


where we have made the change of variable q —)■ q/Tc- It remains to compute the inverse 
Laplace transform on the right-hand side of (25). One hnds three different regimes depending 
on the value of /i. 

1 . 1 < /I < 2 

For 1 < jj, <2 the behavior of s) near s = 1 is given by, (see Eq. (63) in [5]), 


Pi.9, s) - p{9, 1)-(s ^ 1), 

with = l/[asin( 7 r/yu)]. Injecting (26) into (25) yields 

P(a, t) ~ f qiif exp (—) dq (t -q +oo). 

2i7Tr^ Je \rcj 


( 26 ) 


(27) 


Performing then the inverse Laplace transform in (27) and using the reflection formula 
F(— 2 ;)F (;2 + 1) = —tt/ sin( 7 r; 2 ) one obtains, for all 1 < /i < 2, 


sm(7r//i) TTOTc (t/rc)^+'''/^ 

Note that the exponent of t in (28) depends on jj, and (28) tends to (24) as 7 —)■ 1. 


(28) 


2 . 0 < /I < 1 

For 0 < < 1 the behavior of p{g, s) near s = 1 reads, (see Eqs. (63) and (C9) in [5]), 

P{9,s) -PiaA) -^ 2 (^)^i(^)(l-s) (s^l), (29) 


with 


Ma) = 


^+00 


Wi{x) f{g + x) dx, 


(30) 


10 



where wi{x) is defined by its Laplace transform 


^ -A.. A0(A,1) F{k) ,, 

/ Wi[x)e dx = - / dk, 

lo TT Jo 


with F{k) = f{k)/{l — f{k)). Equations (25) and (29) now yield 

^ h{g)Ji{g) f ^ f(lt\ ^ u . . ^ 

p{g,t) --- / exp — dq [t +cx)), 

ZtTT Tr Jr V Tr } 


(31) 


(32) 


which replaces Eq. (28). Performing the inverse Laplace transform in (32) and using the 
reflection formula r(—; 2 )r( 2 : + 1) = — tt/ sin( 7 r 2 :) one obtains, for all 0 < /i < 1, 


sin(7r7)r(l + 7) /2(g)Ji(g) u . . \ 

Pi9,t) -—- (t^+oo). 


(33) 


TTTc (t/rc)l+^ 

Note that the exponent of t in (33) does not depend on p, unlike its counterpart in Eq. (28). 
It follows in particular that (33) does not tend to (24) as 7 —)■ 1. This exponent ‘freezing’ 
from 1 + 7 /p to 1 + 7 as /i decreases past the critical value Pc = 1 (at fixed 0 < 7 < 1 ) is 
due to a switch of leading terms in the asymptotic expansion of p{g, s) — p{g, 1 ) near s = 1 
from cx (1 — s)^/^ for 1 < p < 2, to cx (1 — s) for 0 < p < 1. (See Eqs. (63) and (C9) in [5]). 
Although a term cx (1 — s), as regular, does not contribute to the large / behavior of p{g, /), 
it gives a term oc q'^ in the asymptotic expansion of p{g, 1 — q'^) near g = 0 in (25) which is 
singular if 0 < 7 < 1 and thus contributes to (and dominates) the large t behavior of p(g, t) 
for 0 < p < 1 . 


3. p = 1 and summary 


Finally, at the critical value p = pc = 1, the equation (C9) in [5] gives 

(s -)■ 1 ), 


p{g,s)-p{g,l) - ^^^(l-s)ln^ ^ 


Tra 


1 — s 


which, together with Eq. (25), yields 


P{g^i) ~ [ g^ln(g) exp (— ] dq (t ^ +c 5 o). 


2i7r^aT, 


c j c 


Using then 


1 f y, / ^ f J 1 r (- 7 ) 1 

- / g^ ln(g) exp(gt) dq = — —-r , , ^ ^ 

2m Jc, ^ ’ ^ ’ r(-7) r(-7)2U+7’ 

and the reflection formula r(— 2 :)r (;2 + 1 ) = — tt/ sin( 7 r 2 :), one obtains 

7 sin( 7 r 7 )r(l+ 7 ) / 2 (g) 2 ln(t/re) u . . \ 

Pi9, t) - -7-- - {t ^ + 00 ). 


Tl'^aTr 


(i/r,)‘+T 


(34) 


(35) 


(36) 
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Defining 


(37) 


sin(7r7//i) 

r(i + 7//^), 


and 


Bii = 


7 rsin( 7 r//i) 
sin(7r7) 


TT^ 


r(i + 7), 


we can gather the equations (24), (28), (33), and (36) in a more concise form 




1 

— X 
arc 


Bi 1 - Int( 7 ) < fx < 2, 

< 7^ii72(fi')^(rc/t)^+'^ln(t/rc), /i = 1, 

7rai3ii/2(fi')Ji(fi')(rc/t)^+'^, 0 < /i < 1, 

v 


(t — y H-cxd 
at fixed g) 


(38) 


(39) 


which summarizes the large t behavior of p{g,t) at hxed g and 0 < 7 < 1. The term Int( 7 ) 
in the hrst line of (39) denotes the integer part of 7 and the last two lines are for 0 < 7 < 1 
only. 


IV. ASYMPTOTIC BEHAVIOR OF p(g, t) FOR A LEVY FLIGHT AT LARGE t 
THEN LARGE g 


The random walk is a Levy flight if 0 < /i < 2, in which case the variance of /(?]) does not 
exist and /(t^) ~ at large 7 . In this section we investigate the asymptotic behavior of 

t) for a Levy flight in the limit t —)■ +cx) then g —)■ +cx). 

In the cases 7 = 1 with 0 < /i < 2 and 0 < 7 < 1 with 1 < /x < 2, this asymptotic 
behavior follows readily from Eqs. (24), (28), (36), and the large g behavior of l 2 {g)'- 


h(3) 


^/i/2 


r(l - g/ 2 ) (7M/2 


{g + 00 ), 


(40) 


(see Eq. ( 48 ) in [ 5 ]). One hnds, 

, , 1 Bi fa\^ 

~ r; r(l-^/2)^ (0 (t) (t ^ +00 then g ^ +oo), (41) 

for 7 = 1 with 0 < fj, < 2 and 0 < 7 < 1 with 1 < fj, < 2, and 

P{g,t) -In ^ (t+CX) then 5 (+CX)), ( 42 ) 

arc TT \gj V t J \Tc) 

for 0 < 7 < 1 with /i = 1. 

It remains the case 0 < 7 < 1 with 0 < /x < 1 . Injecting the asymptotics 


0(A,1) ~ 


1 

(aA)^/2 


(A —)■ 0), 


(43) 
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where the integral over x in the second line is equal to r(l — /i/2)r(3/r/2)/r(/i + 1). Thus, 
from (33), (40), and (48) one obtains 

P{9,i) -L +00 then 5 (+CX)). (49) 

arc 2 \g/ \ t / 

Finally, writing Cj = f3ir(l — /r/2)“^, Cu = f 3 ii 7 r“^, and Cm = f 3 iitan(/r 7 r/ 2 )/ 2 , one has 

(a/fi')^(T'cA)^+^/'^, 1 - Int( 7 ) < fj,<2, 

7 C 11 (a/g) {Tc/ty+'^\n{t/Tc), p = I, 

Cm {a/g) (rc/t)^+^, 0 < /i < 1 , 

where Int( 7 ) in the hrst line denotes the integer part of 7 and the last two lines are for 
0 < 7 < 1 only. 
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V. ASYMPTOTIC BEHAVIOR OF p{g, t) FOR LARGE g 

Now, we consider the asymptotic behavior of p{g,t) in the large g limit. If the support 
of the jump distribution f{r]) is bounded, the hrst gap cannot be larger than the (hnite) 
diameter of this support and one trivially has p{g, t) = 0 for every t and g > go with 
go < +c>o large enough. If the support of /{rj) is not bounded the result depends on the tail 
of f{f])- We have singled out three classes of tails which encompass a wide range of jumps 
of practical interest. 


A. Slow decreasing f{r]) 


Slow decreasing jump distributions are dehned by 

f(s + x) ~ f(g) (g ^ +oo), 


(51) 


for any hxed x. Jump distributions with an algebraic tail are slow decreasing. The large g 
behavior of Ii{g,s) is then readily obtained from (17), (18), and (51). One gets 

sf(a) 

h{g,s)r^sf{g) u{y,s)dy = sf{g)(j){0,s) = - (^ ^+c5o), (52) 

h Vl-s 

where we have used 0(0, s) = 1/^/1 — s as given by Eq. (19) in the limit A —)■ 0, (see also 
e.g. [1]). The large g behavior of l 2 {g) depends on f{p) and cannot be specihed further on 
at this point. Thus, Eqs (16) and (52) yield 


^ sf{g)l2{g) ^ ^ 

s)- /-— [g +oo), 


\/l — s 


(53) 


and from Eqs. (20) and (53) one hnds 


p{g,t) = Pia,(g)) exp{qt)dq 


fi.9)h{,g) 

2m 


T(g) 


exp(gt) dq {g —)■ +cx)). 


(54) 


The large t behavior of (54) is then obtained from the small q behavior of '^{q), Eq. (3). 


One has 


1 

2m 


T(g) 


exp{qt) dq 


^\l-^{q) 


2m Jc Q 


exp(qt) , , , 

—dq (t +CX)) 


r(7/2)rc (t/rc)i-T'/2’ 


(55) 
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which together with Eq. (54) gives 

1 f{9)h{g) 




{g —)■ +CX) then t —)■ +cx)). 


(56) 


r ( 7 / 2 ) T , ((/ t ,)>-»/2 

For a Levy flight of index 0 < /i < 2, the large g behaviors of f{g) and h^g) are 
respectively given by (47) and (40), and Eqs. (54) and (56) read 


p{g,i) 


Civr(7/2) 




'At) {g +oo), 


and 

p{g,t) 

respectively, with 

and 


Civ 1-^/2 


arc \g 


Civ = sin(/i7r/2) 


l(t) = - 

^ ^ 2i7r 


{g —)■ +CX) then t —)■ +cx)), 

r(/i + 1 ) 


7rr(7/2)r(l-/i/2)’ 

^(g) 


1 - ^(g) 


exp(gt) dq. 


(57) 

(58) 

(59) 

(60) 


B. Exponentially decreasing /(ry) 

Exponentially decreasing jump distributions are dehned by 

/(g + x) ~/(g)exp(-ca;) (g+cx)), (61) 

for some c > 0 and any hxed x. Substituting (61) for f{g + x) on the right-hand side of (17) 
and using (18), one gets 

hig,s) ^ f{g)s(j){c,s) (g^+cx)), (62) 

and 

hig) ~ f{g) (g ^ +CX)), (63) 

c 

which, together with Eq. (16), yield 

p(g,s) ~ ^^^/(g)^s0(c,s) {g ^+oo), (64) 

and from Eqs. (20) and (64) one hnds 

Pig,^) = exp(gt)dg 

"" ^^(g)0(c,^(g)) exp(gt)dg {g ^+oo). (65) 
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The large t behavior of (65) is obtained from the behavior of (64) near its dominant singu¬ 
larity at s = 1 and the small q behavior of 'h(q'), Eq. (3). The former is given by, (see Eq. 
(60) in [5] with /i = 2 and A = c), 


P{9,s) -P{g,l) - 


0 (c, 1 ) 


n 2 


fidYsy/l — s {g —)■ -|-cx) then s —)■ 1 ). 


( 66 ) 


Injecting then ( 66 ) into the right-hand side of (65) and performing the inverse Laplace 
transform with \l/(g) given by (3), one hnds 




1 U{c,l)yf{gf 

a c 2in 
sin( 7 r 7 / 2 )r(l -|- 7 / 2 ) 

TTOTc 


{Tcqy^'^ exp(gt) dq {g —)■ -|-cx) then t —)■ -|-cx)) 


0 (c, 1 ) 


fig? 


(t/rj 1+7/2’ 


(67) 


where we have used the reflection formula r(—;2)r(2: -f 1 ) = —n/ sin( 7 r 2 :). 


C. Fast decreasing f{r]) 


Fast decreasing jump distributions are dehned by 


f{g + x) ^ f{g)exp{-cg^x + e{x,g)) {g ^+ 00 ), (68) 

for some c, <5 > 0 , any hxed x, and where 9{x,g) is such that \mig^^oo9{x < g~^,g) = 0 . 
Super-exponentially distributed jumps are fast decreasing. Substituting ( 68 ) for f{g + x) on 
the right-hand side of (17) and using (18), one gets, (see Sec. 4.3 in [5] for details). 


h{ 9 ,s) ~ f{g)s<i){cg^,s) 

~ f{g) (s- 


r*+oo 


T^cg^ Jo 


ln[l — sf{k)] dk\ (s' — t +C)o), 


(69) 


and 

( 9 ^+ 00 ), (70) 

eg cg° 

Eqs. (69) and (70), together with Eq. (16), yield 

Pi.g,s) ^ [s - f ln[l - sf{k)]dk'] {g ^+ 00 ), (71) 

cg^ V Jo ) 

from which one readily obtains p{g, 1) for large g, by expanding the logarithm in power series 
of s, (see Eq. (103) in [5]). One hnds in particular, 

P{9, / = 1) ~ {g + 00 ), (72) 

cg^ 


16 



and 


,/) ai 

- rs^ - — f 

, 1) g^oo 

for all I > 1, where ai depends on I only. These results reveal a concentration of p{g, 1) onto 
I = ±1 for fast decreasing jump distributions as g ^ +cx). From this concentration and the 
n +CX) limit of (13), p{g,t) = Eo<«<nP(^IOone has 




P{ 9 ,t) ^ p{g,l = l)p{t\l = 1) {g^+oo), 


(74) 


which reads, using p(t\l = 1) = \h(f) and Eq. (72), 

figfm 


Pi9:t) 


C9° 


{g -P- +oo). 


(75) 


For 7 < 1, T(t) has an algebraic tail at large t, 'F(t) ~ —F(— 7 )r 2 't ^ which follows from 
the small q behavior (3) of T(g), and one gets the large t behavior of (75) as 

sin( 7 r 7 )F(l+ 7 ) /{g^ 


P{9,t) 


{g —)■ +CX) then t —)■ + 00 ), 


TTCTc g^it/TcY^'^ 

where we have used the reflection formula F(— 2 ;)F( 2 ; + 1) = —n/ sin( 7 r; 2 ). 

On the other hand, Eq. (28) with p = 2 and hig) given by Eq. (70) lead to 
sin( 7 r 7 / 2 )F(l+ 7 / 2 ) /{gf 


pigY) 


Trac^Tr 


9^Yt/nY+^/^ 

Equations (76) and (77) suggest the scaling form 


P(9, t) 


with 


where 



acg° 


{t —)■ +CX) then g —)■ +cx)). 


{g —)■ +CX) and t —)■ +cxd). 


Pi^-2(1+1/7)^ (y^O) 

, {y +CX)) 


Vi = sin(7r7) 


r(i + 7) 


(76) 


TT 


(77) 


(78) 


(79) 


(80) 


and 

Vii = sin(7r7/2)^^^i^. (81) 

TT 

Note that ~ ^^9^ corresponds to the cross-over from a ‘concentration’ - or ‘one- 

step’ - regime (for {t/TcY^"^ < acgY where the walker get stuck for a long time t at the 
second maximum and then jumps directly to the first maximum, to a ‘many-steps’ regime 
(for {t/TcY^"^ > acgY where she/he travels a long walk of total duration t (with many steps) 
between the second and the hrst maxima. It remains to check whether or not the scaling 
form (78) really exists. 
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1. Uniform expression ofp{g,t) for large g and t 


Rewriting (19) as 


0(A, s) = 0(A, 1) exp- 


A f+°^ln[l + {l-s)F{k)] 


TT 


A;2 + A2 


dk 


(82) 


in the first line of (69), with F{k) = f{k)/{l — f{k)), one hnds after some straightforward 
algebra, 


P(^,'^(g)) 




cq" 

f{ 9 ? 

cg^ 


'• + 00 


^(g) 


Tiacg'^ 


In I 1 + ^ I dk 
k 


{g —)■ +CX) and g —)■ 0 ), 


(83) 


acg" J 

where we have made the change of variable k = {Tcq)'^^‘^k/a in the integral over k in Eq. (82) 
and use the fact that the integral over k is equal to tt. Inverse Laplace transforming (83) 
w.r.t. g and using (81) yields 

f{g? (V, 


P{9,t) 


cg° 


m 


^ii 


{g —)■ +CX) and t —)■ +cx)). (84) 


Tc{acg^)(t/Tcy+^/^, 

Equation (84) gives a uniform expression of p{g,t) for a fast decreasing f{g) when both g 
and t are large. 

For 7 < 1, one has the large t behavior 'L(t) ~ {t /from which it is readily 

seen that (84) reduces to the scaling form (78) with the scaling function 

1 


K{y) = 


— +Vn 


(85) 


^1 + 2/7 y 

which fulhlls the large and small argument behaviors (79), as it should be. 

For 7 = 1, there is no scaling form such as (78) but the uniform expression (84) makes 

it possible to determine the domains in the (t, g) plane corresponding to the ‘concentration’ 

- or ‘one-step’ - regime and to the ‘many-steps’ regime, respectively. Taking for instance 

4/(t) = exp(—f/Tc) and comparing the two terms on the right-hand side of (84), one 

hnds that the ‘one-step’ regime corresponds to the domain 

5 exp(t/rc^ 
acg > 


2v/T(t/re)3/2’ 

and the ‘many-steps’ regime to the complementary domain (i.e. with > replaced with <). 

Figure 1 summarizes the different behaviors of p{g,t) in the plane {F^‘^,g^) for a fast 
decreasing jump distribution and 0 < 7 < 1. The ‘one-step’ (resp. ‘many-steps’) regime is 
on the left (resp. right) of the diagonal. 


18 



O(g^) 

1 


FIG. 1; Schematic representation of the asymptotic behaviors of p{g,t) for fast decreasing jump 
distributions and 0 < 7 < 1 (amplitudes are given in the text). For i 3> 1 and g < 0(1), p{g,t) is 
given by Eq. (28) with p = 2. For both t S> 1 and g ^ 1, one has the scaling form (78) with K{y) 
respectively given by the first line of (79) ifF^'^ <C g^ and by the second line of (79) if g^. 

For t < 0(1) and g ^ 1, p{g,t) is given by Eq. (75). 

VI. SCALING FORM OF THE ASYMPTOTIC BEHAVIOR OF p{g,t) FOR A 
LEVY FLIGHT AT LARGE t AND g 

In the preceding section we have found that for a fast decreasing jump distribution and 
0 < 7 < 1, the asymptotic behavior of p{g,t) at large t and g takes the scaling form (78). 
Here we consider a Levy flight of index 0 < p < 2 corresponding to a slow decreasing jump 
distribution. We will see that, in this case too, p{g,t) at large t and g can be put into a 
scaling form for all 0 < 7 < 1 . 

Consider hrst the case 7 = 1 . By Eq. (23) and the scaling form (87) in [5], p{g, t) depends 
on t through the combination tg~^ only as both t and g are large. Now, for 0 < 7 < 1, 
the comparison of Eqs. (22) and (25) indicates that t should be replaced with F. Thus, for 
0 < 7 < 1 , p{g,t) is expected to depend on t through the combination Fg~^ only as both 
t and g are large. The simplest scaling ansatz for p{g, t) consistent with this argument and 


I 


D) 


D) 


K(y)~y-2(1+1T) 






/2(g) t 


2 .-1-K/2 
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the asymptotic behaviors ( 50 ) and ( 58 ) reads 

1 f aH^\ 

P{ 9 ,t) -- ^A‘.7 (5--)■+CX) and t-)■+CX)), ( 87 ) 

arc \ 9 j \ 9 ^rJ J 

with the following large and small argnment behaviors for the scaling fnnction: 

Cj 1 _ int(7) < /i < 2, 

~ < Ciiy-^-Fyiny, /r = 1, (y ^ +oo), (88) 

V 

and 

F^,^iy)r^CjyyF^-F\ (y ^ 0). ( 89 ) 

The term Int(7) in the first line of (88) denotes the integer part of 7 and the last two lines 
are for 0 < 7 < 1 only. The scaling form ( 87 ) holds in the regime where both g and t are 
large, with Fg~^ hxed. 

Figures 2 , 3 , and 4 summarize the different behaviors of p{g, t) in the plane {V, g^) for a a 
Levy flight of index 0 < y, < 2. They correspond respectively to the cases 1 — Int(7) < y <2 
with 0<7<1, ya = l with 0 < 7 < 1 and 0 < /i < 1 with 0 < 7 < 1. 


VII. LARGE t BEHAVIOR OF THE MARGINAL DISTRIBUTION Ptime(*) 

We can now determine the large t behavior of the marginal distribution ptime(^) of the 
time between the first two maxima. The marginal distribution Pgap{g) of the gap between 
the hrst two maxima does not depend on 7 and is the same as the one already studied in 
Sec. 3 of Ref. [5] in the discrete time RW setting. 

A. 7 = 1 

For 7 = 1, the large t behavior of ptime(^) is readily obtained by integrating (23) over g 
and by using the large / asymptotics ( 112 ) in [5]. One hnds 

f ^i(t/rc)-i-i /0 l<y<2, 

Ptirae{t) ~ — X < Au ln{t/Tc) {t / Tc)~‘^, h = 1 ) +Oo), (90) 

u 

xliii (t/rc)“^ 0 </i < 1, 
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FIG. 2: Schematic representation of the asymptotic behaviors ofp{g,t) for a Levy flight of index 
1 — Int( 7 ) < pL < 2 and 0 < 7 < 1 (amplitudes are given in the text). For i S> 1 and g < 0(1), 
p{g,t) is given by the first line of Eq. (39). For both t ^ 1 and 5 S> 1, one has the scaling form 
(87) with F{y) given by the first line of (88) iff* ^ g^ and by (89) if E <C g^. For t < 0(1) and 
g ^ 1, p{g,t) is given by Eq. (57). 

with 

Ai = (7ra)-^r(l + l//i) h{gfdg, 

(91) 

^iii = yFuAv) dy- 


B. 0 < 7 < 1 


To determine the large t behavior of ptime(^) for 0 < 7 < 1, fix Ai 3> 1, A 2 = 0(1), and 
for t > (Ai/A 2 )Ot' write 

^+00 ^Ai p+00 

Ptirae{t) = / p{g,t)dg= / p{g,t)dg+ / p{g,t)dg+ / p{g,t)dg. (92) 
Jo Jo JAi JA2t^/i^ 

For 1 < /i < 2 and t —)■ + 00 , (92) is dominated by the hrst two integrals in which one 
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0(g) 

1 


FIG. 3: Schematic representation of the asymptotic behaviors ofp{g,t) for a Levy flight of index 
/i = 1 and 0 < 7 < 1 (amplitudes are given in the text). For t ^ 1 and g < 0(1), p{g,t) is given 
by the second line of Eq. (39). For both t^ I and g^ 1, one has the scaling form (87) with F{y) 
given by the second line of (88) if F ^ g and by (89) if F <C g. For t < 0(1) and g ^ 1, p{g,t) 
is given by Eq. (57) with g = 1. 


Ill 


CN 

LO 


O) 


F(y)-y-'-''>'ln(y) 


/2(g)" r'-''ln(t) 


Oif) 


can use the large t expression (28) of p{g,t). One finds 




sin( 7 r 7 /p) 1(1 +7/p) 




sin(7r/p) naTc{t)Tcy~^^/>^ 
sin(7r7/p) 1(1 + 7 /p) 


h{gfdg 


f +00 


sin(7r/p) TiaTc{tlTcY~^^/^ 


h{gfdg {t + 00 ). 


(93) 


For 0 < p < 1 and t —)■ +cxd, (92) is dominated by the contribution of large, 0{F^y, 
values of g, i.e. by the last two integrals in which one can use the large t and g scaling form 
(87) of p{g,t). One gets 


Ptime(^)-/ 

aTc Ja 


+00 / \ 1 +m( 1 + 1 / 7 ) 

(Jj 




g>^Tc 


dg 


/A') 


++ it/Tcy+^ Jq 

7^111 ln(V+) 

++ it/Tcy+'^ 


dg 

(t —y -{“CX)). 


(94) 
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O(g^) 
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FIG. 4: Schematic representation of the asymptotic behaviors ofp{g,t) for a Levy flight of index 
0 < fi < 1 and 0 < 7 < 1 (amplitudes are given in the text). For t ^ 1 and g < 0(1), p{g,t) is 
given by the third line of Eg. (39). For both t^ 1 and g ^ 1, one has the scaling form (87) with 
F{y) given by the third line of (88) if g^ and by (89) if <C g^. For t < 0(1) and g 1, 
p{g,t) is given by Eg. (57). 





F(y)~y^^2-1/K 

III 

CNI 

CO 

1 

g-i-/^(i+W)F(^K/g/i) 

T 

O) 


. . 



i2(9)M9) r'-'' 


1 


oiin 


where we have made the change of variable y = /g^rf! and used the third line of ( 88 ). 

Finally, for /r = 1 and t —>■ +cxd we use the second line of (88) in the second line of (94) 
with /i = 1 , which gives 


Ptime(^) ~ 


2rc (t/rc)^+^ 


(t —y “t-cx)). 


(95) 


To summarize, for 0 < 7 < 1 the large t behavior of the marginal distribution Ptime(^) is 
given by 


Ptime(^) ~ 


1 

— X 
Tc 


< ln^(t/rc) (t/Tc)"^"^, 
ln(t/rc) (t/rc)“^“^, 

\ 


1 < /i < 2, 
/i = 1, 

0 < /i < 1, 


{t —}■ +CX)), 


(96) 
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with 


^i = o-‘C,r(l-f,/2)2/;“/2(j)2d9, 

A{, = 7"Cr,/2, (97) 

“^ni = 7C111/h) 

where Ci = i3ir(l —/i/2)“^, Cu = jlSu7r~^, and Cm = -Biitan(/i7r/2)/2, the amplitudes Si and 
Su being respectively given by (37) and (38). 

It can be seen that the smaller 7 is, the slower ptime(^) decreases as t increases. The 
third lines of Eqs. (90) and (96) reveal an interesting freezing, as a function of /i, of the 
large t behavior of ptime(^) as fx decreases past the critical value /Xc = 1- It is the continuous 
time counterpart of the freezing phenomenon pointed out in Sec. 5 of [5] for discrete time 
random walks. In particular, it follows immediately from (90) and (96) that the hrst moment 
of Ptiraeit) is uever dehned. This means that although the typical size of t is 0(1), its average 
always diverges with the total duration of the walk. Writing T the total duration of the 
walk, Eqs. (90) and (96) yield the estimate 


1 < /i < 2 


(|t|) ~ p = 1 


T^-'^lnT 0 </i<l 


(T —)■ +CX)), 


(98) 


for 0 < 7 < 1 . 


VIII. NUMERICAL SIMULATIONS 

In this section, we check our main analytical results by comparing them with the ones 
of numerical simulations of the corresponding CTRW. In the simulations, the random jump 
variables piS and the random time intervals between two successive jumps r^’s were drawn 
from the following distributions. 

Jump variables. For p = 2, the jumps were drawn either from a Gaussian distribution 
f{p) = exp(— 77 ^/ 2 )/v^^ - when we investigate the joint PDF of g and t for fast decreasing 
jump distribution - or from a symmetric exponential distribution f{p) = (&/ 2 ) exp (—&|? 7 |) 
- when we study the marginal distribution Ptime(^) “ in which case the amplitudes Ai and 
A!I in ( 11 ) and ( 12 ) can be evaluated exactly. On the other hand, for 0 < /i < 2 , we chose 
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for f{r]) a symmetric Pareto distribution of index fi of the form 


fiv) = 


0 , for |? 7 | < 1 
^ , for |r/| > 1 , 


(99) 


for which it is straightforward to compute the parameter a in Eq. (2). One hnds, 

, 11 //^ 


a = 


r(i-/j)cos(^) 


( 100 ) 


The reason for such a choice is twofold: (i) random numbers drawn from (99) are easy and 
fast to generate numerically (see for instance Ref. [17]) and (ii) f{r]) is a power law for all 
\ri\ > 1 (instead of I?]! 3> 1 ) which allows us to access the asymptotic regime characterizing 
Levy flights rather quickly in numerical simulations. 

Time intervals. Similarly, to generate a CTRW numerically, it is also convenient to draw 
the time variables Tj’s from a distribution T(r) which is a (one-sided) Pareto distribution of 
index 7 < 1 of the form 

{ 0 , for r < 1 

( 101 ) 

7 r ^ , for r > 1 . 

For such a distribution the characteristic time Tc is given by 


r,= [r(l-7)]'/^. 


( 102 ) 


A. Scaling form of the asymptotic behavior of p{g,t) at large t and g 

Computing a joint PDF numerically is notoriously a quite hard task. So, instead of p{g, t) 
we consider the cumulative distribution 

/ OO 

p{g,t’)dt\ (103) 

which yields better statistics at a lesser cost. 

We hrst present the result for a fast decreasing jump distribution f{p) in the sense of 
Sec. VC, i.e. f{p -|-x) ~ f{p)exp{—cp^x) when 7 —)■ cx) [see Eq. ( 68 )]. In this case, one 
expects from Eq. ( 6 ) that in the limit t —)■ ex and g ^ 00 keeping hxed, p>{g,t) 

takes the following scaling form; 

P>(9.*)~a/(9)"(^) ■ (104) 
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where the scaling function H{z) is related to the scaling function K{z) through the relation 

= K(y'l^)dy. ( 105 ) 

From the asymptotic behaviors of K{z) in Eq. (7), one deduces the asymptotic behaviors 
of H{z) as follows 

[ ^ 0 ) 

//(.)- ^ 0 ' ( 106 ) 

I - Z ^ [z ^ CX)) . 

As discussed above, the limit z ^ oo corresponds to the ‘single-step’ regime while the limit 
z ^ 0 corresponds to the ‘many-steps’ regime. 

In Fig. 5 we show a plot of {t /Tc)^'^^'^py{g)[f (g)] ^ as a function of (a c)^/^(t/rc) A( 2 < 5 )^ fQj, 
a fast decreasing jump distribution, namely a Gaussian jump distribution, corresponding to 
(5 = 1, and with 7 = 0.9, for different values of t = 5000,10000, 20000 and 40000, while the 
total duration of the CTRW is T = 10^. The data have been obtained by averaging over 
10® independent realizations of the CTRW. We see on this plot that, for sufficiently large 
values of g, the curves for different times t tend to collapse on a master curve, as predicted 
by Eq. (104), which indeed holds only in the limit g ^ 00, t ^ 00 with hxed. For 

smaller values of g we see a clear deviation from this scaling form. We emphasize on the fact 
that, in our simulations, we have been able to access the regime where the scaling variable 
CX is small (corresponding to the ‘many-steps’ regime), where H{z) is expected 

to behave as H(z) cx z~‘^. It is extremely difficult to reach the complementary ‘single- 
step’ regime corresponding to 3> 1. Indeed, for such high values of t > 5000 the 

condition ^ ^ with (5 = 1 and 7 = 0.9 yields g 3> 46, which is already large enough 

to prevent any access to reliable statistics (to understand the difficulty, just realize that 
for the Gaussian jump distribution we consider in the ‘single-step’ regime, the typical gaps 
are greater than 10 less often than one time out of the Avogadro’s number). Gonsequently, 
the clear numerical observation of the crossover behavior between the ‘many-steps’ and the 
‘single-step’ regimes described by the scaling function in Eq. (104) goes beyond the scope 
of the present paper. Notice that, for 7 = 1 and t ~ Tc, the ‘single-step regime’ was clearly 
demonstrated in Ref. [5]. 

We now present numerical results for Levy flights (0 < p < 2), for which /(?]) is slow 
decreasing in the sense of Sec. V A. In this case, one expects from Eq. (8) that in the limit 


26 




FIG. 5: Plot of {t/Tc)^'^^‘^p>{g)[f{g)]‘^ as a function of {ac)^^^{t/Tc)~'^^^‘^^'^g for a fast decreasing 
jump distribution, namely a Gaussian jump distribution, corresponding to <5 = 1, and with 7 = 0.9. 
The different curves correspond to different values of t = 5000,10000, 20000 and 40000 for a CTRW 
of duration T = 10®. The (relatively) good collapse of the different curves for large values of g 
corroborate the scaling form predicted in Eq. (104). The solid line is a guide to eyes, corresponding 
to the small argument behavior of the scaling function H{z) in Eq. (106) with <5 = 1. Inset: same 
data as in the main panel but without the rescaling. 


t —>■ 00 and —)■ 00 at fixed gt t) takes the following scaling form: 
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where the function is related to the scaling function through the relation 

1 r 




yl+/^ 


^1/W^ 
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The small and large behaviors of G^^^{z) are readily obtained from the asymptotic behav¬ 
iors of Ffj,^^{z) given in Eqs. (9) and (10). One gets 
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for large z, where C'm = dy, and 


GiJ.,-y{z) ~ < 


7 

^iii z ^ , 0 < /i < 1 . 

for small z, where the last two lines are for 0 < 7 < 1 only. 


Ciz ^ , 1 — Int(7) < fi <2 
—ln(l/2;) , /r = 1 (;2 —>■ 0) , 


( 110 ) 



FIG. 6 : Plot of a P>{g^ t) as a function of {g/a){Tc/t)'^^^ with 7 = // = 1/2 for different 

values of t = 10,100 and t = 1000, for a CTRW of total duration T = 10^. The good collapse of 
the curves for different values of t confirm the scaling form predicted in Eq. (107). The dashed 
and solid lines are guides to the eyes, indicating the algebraic behaviors predicted in Eqs. (109) 
and (110). Inset: same data as in the main panel but without the rescaling. 

In Fig. 6 we show a scaled plot of a p>(( 7 , t) as a fnnction of {g/a){Tc/t)"'^^ 

with 7 = /i = 1/2 for different valnes of t = 10,100 and t = 1000, for CTRWs of total dnra- 
tion T = 10^. The data have been obtained by averaging over 10® independent realizations 
of the CTRW. This plot demonstrates a good agreement with onr predictions for the scaling 
form of the joint PDF in Eq. (107) as well as with the asymptotic behaviors predicted in 
Eqs. (109) and (110). 

28 













B. Asymptotic behavior of ptime(t) at large t 


An interesting feature of our results for the distribution Ptime(^) of the time t between the 
hrst two maxima is the presence of logarithmic corrections for 0 < /i < 1 and 0 < 7 < 1 
[see Eq. (12)]. To demonstrate these logarithmic corrections numerically, it is convenient to 
study the distribution piogtime(E) of the logarithm of the time, Y = ln(t/rc) instead of the 
time distribution itself. From piogtime(E) = Tce^ptime('rce^) and Eq. (12) it follows that the 
asymptotic behavior of piogtime(l^) at large Y and for 0 < 7 < 1 is given by 


Plogtime(^) ^ \ 


A!i exp , 1 < /i < 2 

A'nY^ exp(- 7 F),p = l 
A'ljj Y exp (—7 F) , 0 < /i < 1 , 


with the amplitudes 


A = (77r/2)r(l + 7/2) 

= ^r(l + 7)sin(77r) , 


( 111 ) 


( 112 ) 

(113) 

(114) 


where the value of Aj corresponds to a symmetric exponential jump distribution (see the 
beginning of this section). 

In Fig. 7, we show a plot of piogtime(E) as a function of F for three different values of p; 
p = 2 (with exponentially distributed jumps), p = 1 , and p = 1 / 2 , that are representative 
of the three different regimes in Eq. (Ill), for a hxed value of 7 = 1/2. This plot shows that 
the large F behavior of piogtime(l^) does depend on p as predicted in Eq. (111). To make 
a precise comparison between our numerical results and our analytical predictions (and in 
particular to demonstrate the expected logarithmic corrections), we now analyze these data, 
for different values of p, more carefully. 

Fig. 8 corresponds to the case of a sub-diffusive walk with p = 2 (with exponentially 
distributed jumps) and 7 = 1/2. We show a plot of piogtime(l^) exp( 7 F//i) for 7 = 1/2 as 
a function of F = log(t/rc) for CTRWs of increasing length, namely T = 10^, T = 410^ 
and T = 10®. Here, the data have again been obtained by averaging over 10® independent 
CTRWs. It can be seen that in the time interval 1 t T where our prediction actually 
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FIG. 7: Plot of piogtime(^) as a function oi Y = ln(t/rc) on a log-linear plot for 7 = 1/2 and for 
different values of /r = 2,1 and 1/2, corresponding to the three different regimes in Eq. (111). The 
solid lines are guides to the eyes corresponding to the asymptotic behaviors for ^ = 2 and n = If2 
given in Eq. (111). These data have been obtained for CTRWs of length T = 10®. 

holds, piogtime(^) exp( 7 F//i) approaches a constant which gets closer to the expected value 

as the total duration T of the walk increases. The observed small upward discrepancy, as 
well as the steep upward bending for larger Y, are hnite size effects. Note the very stretched 
vertical scale. For the longest walks with T = 10®, numerical and analytical results are in 
agreement to within an accuracy of less than 5%. 

Next we consider the case /x = 1/2, which corresponds to the third regime in Eq. (111). 
In Fig. 9 we show a plot of piogtime(F') exp( 7 ’F), for p = 1/2 and 7 = 1/2 as a function 
of Y and for different durations of the CTRW, T = 10^, T = 4 ■ 10^, and T = 10®. As 
before, the data have been obtained by averaging over 10® independent CTRWs. The linear 
behavior, in a good agreement with our analytical prediction oc R, is a clear signature 
of the logarithmic correction of ptime(t) as expected from our analytical result in Eq. (12). 

Finally, Fig. 10 shows a plot of piogtime(R) exp (qR) as a function of R, for p = 1, 
7 = 1/2, and the same walk durations as before (averaging procedure is also the same). It 
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FIG. 8; Plot of piogtime(^) exp( 7 y/|u) for 7 = 1/2 and /r = 2 (namely exponentially distributed 
jumps) as a function of y = log(t/rc) for CTRWs of increasing length, namely T = 10^, T = 4 • 10^ 
and T = 10®. Note that the upward bending is a finite size effect, as our predictions actually hold 
for 1 <C t <C T which is confirmed by the fact that this bending happens for larger and larger 
values of t as the total duration T of the CTRW increases. In the appropriate regime such that 
1 t <C T Piogtime(^) exp( 7 y///) Converges to a constant which is indeed close (with an accuracy 
of less than 5% for T = 10®) to our analytical predictions in Eq. (111). 

can be seen that onr nnmerical data are compatible with the expected qnadratic behavior 
(111), Plogtime(h') exp (yy) r\j AliiY"^. Recalling that Y = log(t/rc), this qnadratic behavior 
conhrms onr analytical predictions for Ptime(^) in Eq. (12). 


IX. SUMMARY AND PERSPECTIVES 

In this paper, we have performed a detailed analytical stndy of the statistics of the gap 
g and the time interval t between the hrst two maxima of a CTRW in the limit where 
the nnmber of steps in the walk goes to inhnity. The resnlts we have obtained are qnite 
general as we have addressed the qnestion for the wide class of a symmetric, bonnded, and 
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FIG. 9: Plot of piogtime(^) exp ( 7 I") as a function of Y for /r = 1/2 and 7 = 1/2 and for different 
durations of the CTRW: T = 10*^, T = 4 • 10^, and T = 10^. The solid line, which corresponds 
to our analytical result in Eq. (Ill) with = 0.0448968... is a guide to eye, indicating a good 
agreement between our analytical results and numerics. This linear behavior is a clear signature of 
the logarithmic corrections inptiine(i)) which is a peculiarity of the CTRW model. Upward bending 
at large T is a finite size effect which decreases as T increases. 

(piecewise) continuous jump distribution, including the case of Levy flights of index /i, with 
0 < yU < 2. The average trapping time (r) between two successive jumps can be either 
hnite, (r) = Tc < +cxd, or inhnite. In the latter case, the trapping time distribution T(r) is 
assumed to have a power law tail ^(r) oc with 0 < 7 < 1. 

For a hnite average trapping time between two successive jumps, our results coincide 
with the ones in [5], devoted to discrete time RWs, in which the number of jumps between 
the hrst two maxima, I, is merely replaced with t/Tc, as expected from simple law of large 
number arguments. On the other hand, our study provides a non trivial extension of the 
one in [5] to the case where the average time between two successive jumps does not exist 
and law of large number arguments cannot be used. In particular, by taking 7 < yu/2 it 
makes it possible to address the question for sub-diffusive RWs, out of reach of the discrete 
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FIG. 10: Plot of piogtime(^) exp ( 7 y) as a function of F for /i = 1 and 7 = 1/2 and for different 
durations of the CTRW, T = 10^, T = 4 • 10^, and T = 10®. The solid line, which corresponds to 
our analytical result in Eq. (Ill) with Ain = 0.000357277... is a guide to eye, indicating a good 
agreement between our analytical results and numerics. Upward bending at large U is a finite size 
effect which decreases as T increases. 

time analysis of [5]. 

We have first shown that the joint PDF Pn{g, t) converges to a stationary PDF, Pn{g, t) —)■ 
p{g, t), as n —)• oo, where n is the total number of steps in the walk. We have then obtained 
an explicit expression of the Laplace transform of p{g,t) w.r.t. t, given in equations (20) 
and (17), from which we have performed a detailed analysis of p{g,t) in the plane {g,t) 
p and 7 in the whole ranges 0 < p < 2 and 0 < 7 < 1 , and for the three different main 
classes identihed in [5]: (i) slow, (ii) exponentially, and (hi) fast decreasing f{p) at large p. 

In the case (hi) of a fast decreasing jump distribution (/i = 2) with 0 < 7 < 1 and in 
the scaling regime g, t ^ 1 with hxed we have shown that p{g,t) takes a scaling 

form, Eqs. ( 6 ) and (7). The switch from the first to the second behavior (7) corresponds to 
the cross-over from a ‘concentration’ - or ‘one-step’ - regime where the walker get stuck for 
a long time t at the second maximum and then jumps directly to the hrst maximum, to a 
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‘many-steps’ regime where she/he travels a long walk of total duration t (with many steps) 
between the second and the hrst maxima. We have summarized these results in Figure 1. 
The numerical observation of the cross-over and ‘one-step’ regime is still out of reach as it 
requires a prohibitively large value of g. Note, however, that we expect a similar scaling form 
for p{g,t) if f{ri) has a bounded support —r/max < ''7 < Vraax and g —)■ -|-cx) is replaced with 
g —^ ^max, the largest possible (bounded) value of the gap. This should make a numerical 
observation of the cross-over and ‘one-step’ regime easier in this case. The study of p{g, t) for 
jump distributions with a bounded support (in both discrete and continuous time settings) 
will be the subject of a future work. 

For type (i) Levy flights of index p, with 0 < p < 2, and in the scaling regime g, t ^ 1 
with hxed t'^g~^, we have shown that p{g, t) takes the scaling form ( 8 ) with the asymptotic 
behaviors (9) and (10). These results, summarized in Figures 2 to 4, provide a non trivial 
extension of the ones in [5] to the case 0 < 7 < 1. 

Finally, from the stationary joint PDF p{g,t) we have computed the stationary marginal 
distribution, Ptime{t), of the time between the first to maxima of the walk. Its asymptotic 
behavior, which depends on 7 and the Levy index p, is summarized in equations (11), (12), 
and in Table 1. The third lines of (11) and (12) reveal that the freezing phenomenon of 
the large t behavior of ptime(^), as a function of p, as p decreases past the value /ic = 1 , 
which was first put in evidence for discrete time RWs in [4, 5], also exists for CTRWs. As 
a consequence, the first moment of Ptime(^) is never dehned. 

Our study raises interesting open questions. For instance, this work focuses on the hrst 
gap and it would be interesting to study higher order gaps. In Ref. [11], the study of the k-th. 
gap of long random walks, for p = 2 revealed a rich and universal behavior (i.e., independent 
of the jump distribution) of the stationary PDF of the k-th gap in the large k limit. The 
extension of these results, which also hold for CTRW with p = 2 and any value of 0 < 7 < 1, 
to Levy hights, i.e. 0 < /i < 2, remains a challenging open question. Similarly, it would 
very interesting to study the time between the k-th and {k l)-th maximum for CTRW 
(both for 7 = 1 and 7 < 1), for which one also expects universality for /c 3> 1. Finally, the 
generalization of our results to higher dimensions would also be very interesting. 


34 



[1] Comtet, A., Majumdar, S. N.: Precise asymptotics for a random walker’s maximum. J. Stat. 
Mech. P06013 (2005). 

[2] Derrida, B., Brunet E.: Statistics at the tip of a branching random walk and the delay of 
traveling waves, EPL 87 60010 (2009). 

[3] Derrida, B., Brunet E.: A branching random walk seen from the tip, J. Stat. Phys. 143 , 420 

( 2011 ). 

[4] Majumdar, S. N., Mounaix, Ph., Schehr, G.: Exact Statistics of the Gap and Time Interval 
Between the Eirst Two Maxima of Random Walks, Phys. Rev. Lett. Ill, 070601 (2013). 

[5] Majumdar, S. N., Mounaix, Ph., Schehr, G.; On the gap and time interval between the first 
two maxima of long random walks. J. Stat. Mech. P09013 (2014). 

[6] Schehr, G., Le Doussal, P.; Extreme value statistics from the Real Space Renormalization 
Group: Brownian Motion, Bessel Processes and Gontinuous Time Random Walks, J. Stat. 
Mech. P01009 (2010). 

[7] Bouchaud, J.-P., Mezard M.: Universality classes for extreme-value statistics, J. Phys. A 30 , 
7997 (1997). 

[8] Dean D. S. , Majumdar S. N.; Extreme Value Statistics of Hierarchically Correlated Variables: 
Deviation from Gumbel Statistics and Anomalous Persistence, Phys. Rev. E 64, 046121 (2001). 

[9] Pollaczek E.: Eonctions caracteristiques de certaines repartitions definies au moyen de la 
notion d’ordre, C. R. Acad. Sci. Paris 234, 2334 (1952). 

[10] Spitzer P.: A combinatorial lemma and its application to probability theory, Trans. Am. Math. 
Soc. 82 , 323 (1956). 

[11] Schehr, G., Majumdar, S. N., Universal Order Statistics of Random Walks, Phys. Rev. Lett. 
108 , 040601 (2012). 

[12] Scher, H., Montroll, E. W., Anomalous transit-time dispersion in amorphous solids, Phys. 
Rev. B 12, 2455 (1975). 

[13] Bouchaud J. P., Georges A., Anomalous diffusion in disordered media: Statistical mechanisms, 
models and physical applications, Phys. Rep. 195 , 127 (1990). 

[14] Metzler, K., Klafter, L, The random walk’s guide to anomalous diffusion: a fractional dynamics 
approach, Phys. Rep. 339 , 1 (2000). 


35 



[15] Carmi, S., Turgeman, L., Barkai, E., On distributions of functionals of anomalous diffusion 
paths, J. Stat. Phys. 141, 1071 (2010). 

[16] Franke, J., Majumdar, S. N., Survival probability of an immobile target surrounded by mobile 
traps, J. Stat. Mech. P05024 (2012). 

[17] De Mulatier C., Rosso A., Schehr G., Asymmetric Levy flights in the presence of absorbing 
boundaries J. Stat. Mech. P10006 (2013). 


36 



